Dynamically Error-Corrected Gates for Universal Quantum Computation 
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Scalable quantum computation in realistic devices requires that precise control can be implemented 
efficiently in the presence of decoherence and operational errors. We propose a general constructive 
procedure for designing robust unitary gates on an open quantum system without encoding or 
measurement overhead. Our results allow for a low-level error correction strategy solely based on 
Hamiltonian engineering using realistic bounded-strength controls and may substantially reduce 
implementation requirements for fault-tolerant quantum computing architectures. 
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Physical realizations of quantum information process- 
ing hold the potential to solve problems in physics simula- 
tion, combinatorial analysis, and secure communications 
virith unprecedented power compared to known classical 
counterparts []]. The discovery that arbitrarily accurate 
(fault-tolerant) quantum computation (QC) may be sup- 
ported by real-world imperfect devices provided that the 
overall noise is below a certain threshold [2] indicates 
that no fundamental obstacle prevents this power to be 
harnessed in principle. The requirement on errors affect- 
ing an operation on physical qubits may be quantified in 
terms of an appropriate error per gate, EPG. While sug- 
gestive evidence exists that EPGs well above 10~^ can be 
tolerated through exploitation of concatenated quantum 
error correction (QEC) [>>] and post-selection [4, ^>], im- 
plementation requirements remain daunting. The main 
obstacle in generating precise unitary transformations on 
an open quantum system is due to the fact that a control 
prescription realizing a desired gate in the ideal (closed- 
system) limit no longer works accurately in the presence 
of an uncontrollable environment: typically, the resulting 
EPG will be proportional to both the noise strength and 
gating time. Our goal is to show how, for arbitrary open 
quantum systems undergoing linear decoherence, infor- 
mation about the errors of a fixed set of primitive gates 
may be exploited to construct dynamically corrected gates 
(DCGs) which achieve a significantly smaller EPG using 
realistic bounded-strength control resources. 

Our approach successfully merges elements from differ- 
ent techniques for high-fidelity coherent quantum control 
- including composite [G, 7, 8] and strongly-modulating 
pulses [9] from nuclear magnetic resonance (NMR), as 
well as dynamical decoupling (DD) methods for decoher- 
ence suppression [10, 11]. Common to these approaches 
is the idea of enforcing active error cancellation through 
purely Hamiltonian (open-loop) control, bypassing the 
need for measurement or feedback implicit in QEC. Sim- 
ilar to strongly- modulating pulses, DCGs coherently av- 
erage out unwanted evolution by cascading primitive con- 
trol operations within a self-contained composite block. 
Unlike the standard NMR setting, however, where the 
unintended error component is either classical or induced 



from a known spin Hamiltonian, DCGs must operate 
without assuming complete knowledge or control over the 
underlying open-system Hamiltonian. As we shall see, 
a general analytic prescription for the required control 
modulation may be established by suitably incorporating 
known relationships between errors into control design. 

Our results advance open-loop approaches to error con- 
trol in several ways. While DD-inspired constructions 
have been employed to obtain shaped pulses which ap- 
proximate ideal '(5-pulses' by sclf-refocusing specific un- 
wanted couplings up to some degree [12, 13], existing 
schemes do not incorporate the effect of a generic quan- 
tum environment. In contrast, DCGs provide a complete 
prescription for achieving universal QC with reduced er- 
ror solely based on unitary manipulations. Although pro- 
cedures for combining DD with logic gates have been es- 
tablished in [14] in the so-called 'bang-bang' limit and 
rigorously analyzed in [15], our approach has the key ad- 
vantages of avoiding unphysical unbounded controls from 
the outset, along with the need for encoding and stringent 
time-synchronization. Thus, DCGs can be instrumental 
as a low-level error control strategy for bringing quantum 
fault-tolerance closer to current capabilities. 

Error and control assumptions.— Consider a target sys- 
tem S consisting of n qubits, coupled to an environment 
(or bath) B. The joint evolution is governed by a total 
Hamiltonian of the form H = Hs (E) Ie + H^, where the 
error Hamiltonian — Is ® Hb + Hsb is responsible 
for both pure bath evolution via Hb and unwanted in- 
teraction via Hsb = X]q '^a ® Ba, for operators Sa{Ba) 
acting on S (B), with Sa traceless. We focus on the case 
where S is driftless. Hs ~ 0, and subject to arbitrary 
linear decoherence: 
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We further assume the bath operators Hb and Ba ' 
bounded [Ki], but otherwise unknown. 

Control over S is implemented by applying a time- 
dependent Hamiltonian Hctr\{t)- Let the propagator 
Uctr\{t2-,ti) = 7I|.exp{-i J^^ Hctri{t)dt}, in units ?i = 1, 
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effect an intended gate on S in tlie absence of H^. If 
HsB 7^ 0; application of the same control Hamiltonian 
results in an actual propagator U{t2,ti) whose action 
deviates from the intended one due to error dynamics 
induced by H^. The EPG may be quantified in terms 
of an Hermitian error phase operator ^(t2,ti) by writ- 
ing U{t2,ti) = C/ctri(i2,<i)exp[-i$(<2,ii)]- Physically, 
$ is related to the time-dependent error Hamiltonian 
which describes the joint system-bath evolution in the 
'toggling frame' that follows the control [10]. The norm 
of $ bounds the fidelity loss between the intended and 
actual evolution of the system [17]. $ may contain pure- 
bath terms that have no effect on the reduced system 
dynamics, e.g., when Hsb = 0, $(t2,ii) = (^2 — ti)HB. 
To avoid ambiguity, we define an operator A modulo pure 
bath terms by letting A mods = A — 2~"TrsA. 

We specify the available control resources by assuming 
access to the following switchable control Hamiltonians: 

{h, {t)x(^ , hy , {t)z^'^ z^-'") }, (2) 

for appropriate control inputs ha (t) . This allows any uni- 
tary evolution on 5* to be efficiently approximated within 
the circuit model of QC [1]. Specifically, a universal 
gate set is given by: (i) NOOP gates, in which no oper- 
ation is performed on some or all of the qubits; (ii) ar- 
bitrary single-qubit rotations on qubit i, X^g'' = g-'9A'<^' 
and Y^p = e~**^'^' ; (iii) two-qubit phase gates between 

qubits i and j, Z^^g'' = g^*^^*'^*^'. Realistic control 
profiles ha{t) will be constrained in many ways due to 
limited pulse-shaping capabilities. We incorporate finite- 
power and finite-bandwidth constraints by assuming the 
existence of a minimum switching time r^in for modu- 
lation and by requiring all control strengths ha{t) to be 
bounded by /imax- A gate realized using a single control 
input in a pre-determined manner shall be referred to as 
primitive - for instance, gates implemented by turning 
on and off a Hamiltonian from the above set according 
to a rectangular profile may be called primitive. 

In control-theoretic terms, our goal is as follows: Given 
a desired unitary gate C/gatc in the universal set, devise 
a control procedure -ffctri(i) for effecting a DCG using 
Hamiltonians in the available repertoire, Eq. (2), such 
that (i) Ucti\{t2,ti) = C/gate; (ii) the error \\^{t2,ti) 

mods 1 1 

is significantly reduced compared to the primitive EPG. 
We construct an analytic perturbative solution which 
does not resort to measurements, extra qubit resources, 
or any quantitative knowledge of the error Hamiltonian 
iJe except its algebraic structure, Eiq. (1). Our solution is 
perturbative in the sense that ||$(t2, ti) mods || becomes 
proportional to r^j^. Since naive switching of -ffctri(0 
would yield error phases that scale with 0(rmin ||i^ss||)j 
this implies EPGs reduced by a factor of 0(Tmin ||^^e||)- 
Error Combination and Cancellation.— The first step 
toward DGCs is to quantify the error phase arising from 
cascading N gates. Let U = UnUn-i ■ ■ ■ UiUo, where 
Uq = IsB and the j-th operation Uj during [<j_i,tj]. 



which is intended to generate C/ctrij = Uctri{tj ,tj-i) in 
the absence of He, has an EPG ^j. Up to the first order 
in maxj(||$j||), the total error $[/ is: 

N 

•^c/ = E uU.^^'^jU.t.ij-i + $['+1 , (3) 

where H^'^+'H = 0(A^^ max ||^). Each error $j can 
be computed up to the first order in {ti — ||i?sB|| as 

$™ = / ' Ul^.it, t,.i)HeU,ui{t, tj.i)dt + $f +1, 

with = 0[(t,-i,_i)2 \\Hef]. Here, Al^+l includes 

corrections of second- or higher-order powers in to A. 

The next step is to seek a combination of gates which 
removes the combined error (at least) up to the first 
order. This is reminiscent of DD approaches to QEC, 
whereby the time scale difference between the action 
of the errors (in the non-Markovian regime) and the 
available controls is leveraged to reduce or symmetrize 
the effect of the environment [10]. While several fla- 
vors of DD arc possible depending on system and de- 
sign specifics, the formulation relevant to our purpose 
is Eulerian DD (EDD) [IS], which implements DD us- 
ing bounded-strength controls and guarantees robustness 
against systematic control faults. Consider a set of uni- 
tary operators on S, which form a (projective) represen- 
tation {Ug-} of the so-called DD group Q ~ {.9i}|Li, and 
let VI be the subspace of all traceless (modulo pure bath) 
operators which obey the following decoupling condition: 

f y t/] =0, \/Een. 

\^ y» «V mods 
i 

A good DD group ensures, in particular, that H contains 
all operators generated by the errors {So} we wish to 
correct. For single-qubit error generators as in Eq. (1), 
the smallest DD group is 5 = Z2 x Z2 under the n-fold 
product representation in terms of collective Pauli matri- 
CCS {/(an)^x(-"),y(-"),Z(-")}, with Sf'^ = ®'US^S'. 
Consider the Cayley graph associated to a set of gener- 
ators {hj}'-^i of Q. In this directed graph, each vertex 
represents a group element, and two vertexes gi, gj are 
connected with an edge labeled by the generator h if and 
only if gj = gih. Cayley graphs always have an Eulerian 
cycle, that is, there exists a closed sequence oi LxD con- 
nected edges that visits each edge exactly once. The Cay- 
ley graph of Z2 X Z2 is shown in Fig. l(left), along with 
an Eulerian cycle beginning (and ending) at the identity. 
Consider now a sequence of primitive gates that imple- 
ment the generators {hj} in the same order they appear 
in this Eulerian cycle. Provided that same generators are 
implemented by same gates along the path, the error of 
the full sequence may be obtained using Eq. (3): 

hi 



^^^^^ .^<r^^v, 



Figure 2: Quantum circuit for a DCG on qubit 1. The NOOP 
implementation is recovered by dropping all C±20 gates. 



Figure 1: Left: The Cay ley graph for Z2 x Z2 represented 
by collective Pauli operators {I , X,Y, Z} , together with an 
Eulerian cycle marked by numbers. Edges are labeled by the 
generators X, Y, and arrows denote the direction of action. 
Right: Modified Cayley graph supporting an Eulerian path 
that cancels errors in a non-identity gate. The new edges 
correspond to sequences with the same (leading-order) error. 



where is the EPG associated to hj and ||$eddII ^ 
0(max ||<i>/i^. IP). As long as G up to the first or- 
der, then irrespective of how each primitive gate is im- 
plemented, ($EDD)modB = up to the fiist Order. Thus, 
a dynamically corrected NOOP can be effected with a sig- 
nificantly smaller error compared to the free evolution. 

Dynamically Corrected Gates.— The above construc- 
tion does not directly extend to transformations other 
than NOOP [lit]. In order to build composite gates which 
(as in EDD) not only cancel the combined error but (un- 
like EDD) effect a non-trivial rotation U, we need to use 
relationships between the errors of the primitive gates. 
One such relationship is obtained when two combinations 
of gates Mi and Mu, intended to generate / and U I, 
respectively, have the same error to the leading order. 
Exphcitly: Mu = U exp{-i^) and Mj = exp(-i<i>). 
This relationship suggests the following modifications of 
the Cayley graph used for NOOP: (i) To every vertex other 
than identity, attach a self-directed edge labeled with Mj; 
(ii) To the identity vertex attach a new vertex represent- 
ing U through an edge labeled by Mu, see Fig. 1 (right). 
This new graph possesses an Eulerian path starting at / 
and ending at U, which implements a DCG U. The net 
EPG for the corresponding gate sequence is given by: 



$DCG = $EDD + Ul^Ug, + $ 



[2+] 
DCG' 



where 11$ 



1-2+] 
DCG I 



0[max(||$^JlM|$||2)]. As long as 
and <I> all belong (up to the first order) to the subspace 
ft of correctable errors, ($DCG)mod_B = 0, up to 0(||$|[). 
Thus, U is implemented with a significantly smaller error 
compared to Mu- 

Resource Requirements and Performance.— We now 
provide an explicit construction of dynamically corrected 
NOOP, single-, and two-qubit gates within the control sce- 
nario specified in Eq. (2). Let 6 parametrize the rotation 
angle. For each operation, we assume that: (i) the control 
profile hg(t) in the interval [ti,t2] is obtained by stretch- 



ing, scaling, and additions of a fixed reversible pulse 

shape ho{t), t e [0, 1], e.g., we let 4*''*'' W = ^hoi^J; 
(ii) both positive and negative values for holt) are avail- 
able. Rectangular pulses provide the simplest illustrative 
setting, with holt) = 1 if i G [ii,i2]7 and zero otherwise. 
A concrete example of control sequences that yield equal 
errors, yet differ in the intended action is given by the 
following two piecewise-constant control profiles: 

hi{t) = 0ho{t/T)-eho[l-{t-T)/T], (4) 

h2{t) = 0ho{t/2T), r>T,„i„. (5) 

In Eq. (4) , hi corresponds to a sequence of two primitive 
gates intended to implement the identity. Notice the re- 
verse of the basic shape as ho[l — {t — t)/t]. In Eq. (5), 
/i2 corresponds to a primitive gate of duration 2t imple- 
menting Ce. One may show [ 1 0] that for any choice of the 
basic shape ho, the errors are the same up to the leading 



order, $ 



[1] 



$2 , and belong to the subspace Q2.i 



of inhomogeneous two-qubit system-bath terms spanned 
by {a^a^a'-J^ ® B^^p}, with a ^ p. The modified Cayley 
graph for a DCG C2e can be obtained by specializing the 
Eulerian path depicted in Fig. 1 to the case where U 
originates from the control sequence /12 and / from hi, 
respectively. Notice that the final edge associated with 
h2{t) connects the identity to the new vertex represent- 
ing the desired C2e- The required gates x'''''' and F^*^") 
arc implemented through collective single-qubit Hamil- 
tonians. Direct calculation shows that the errors for all 
the primitive gates used involve only single-qubit terms, 
^^singio C ri2,inhom C il. Thus, thc combincd error is zero 
modulo pure bath terms, up to thc first order. The re- 
sulting DCG circuit appears in Fig. 2. The DCG for €26 

is longer in duration than a single primitive C2e gate by a 

[2+1 

factor of 16. The error || (<&ocg) mods || on the other hand 
is found to be reduced by a factor of order 0{t \\He\\). 

The above estimate indicates that significant improve- 
ment is expected for sufficiently small Tmin- Quantitative 
supporting results are given in Fig. 3 for an algorithm 
which prepares a 3-qubit cat state {tpcat) for a spin-based 
qubit device coupled to a spin bath. While the simulation 
is performed for a relatively small system, it does show 
the expected behavior of DCGs in the regime where in- 
dividual primitive EPGs are small. For this reason, the 
net fidelity loss using primitive EPG is chosen as the x- 
axis in Fig. 3, as opposed to a less objective measure 
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Figure 3: DCG performance in a 3-qubit cat-state algorith- 
mic benchmark, using the metric 1 — -^pcatPoutPcat, with 
Peat = \ip cat) {ip catl and pout the actual output state. Lin- 
ear decoherence from a 5-qubit bath and systematic control 
errors are included. Hadamard and CNOT gates appearing in 
the inset quantum circuit are decomposed as sequences of (2 
and 6, respectively) physical primitive gates, and each such 
gate is replaced by a DCGs, for a total of (2 -I- 2 * 6) * 16 = 
256 control time slots. Every bath spin la interacts via a 
Heisenberg coupling with every other system or bath spin, 
He=r J2a<b ■P'''^ +A Y,^,a ^^^^ ' ' ^uch a way that 
(in units of t^I^ F = 1 and logjQ ^ — — 1,-1.4, ••• ,—5.8, 
corresponding to different circles along a given curve in the 
figure. Rectangular pulse profiles are used throughout, sys- 
tematic pulse- length errors being included by modifying ho(t) 
as /io(i)(l + e)- Line thickness proportionally represents error 
strength e. Note the bunching of data points at low EPG for 
each e, signaling a regime dominated by systematic error. 



such as the coupUng strength. A systematic control er- 
ror is also included, allowing its interplay with the bath- 
induced errors to be analyzed. Several conclusions may 
be drawn. For ideal control, the EPG reduction is mani- 
fest in the change of slope between dynamically corrected 
(thick dashed line) and uncorrected (narrow dashed line) 
data, resulting in the 'cone of improvement' of the proce- 
dure. The smaller the primitive EPG, larger amounts of 
systematic control faults may be intrinsically tolerated, 
resulting in larger regions of improvement within the 
above cone. Notice that while systematic error compen- 
sation along the standard Cayley graph is ensured [18], 
robustness need not be retained along the added arms in 
the modified graph. Thus, no further improvement from 
reducing the primitive EPG arises once uncompensated 
systematic error dominates over bath-induced error, lead- 
ing to the observed performance plateaux. 

Discussion.— We have shown how to synthesize unitary 
gates able to approximate ideal gates in a universal set 
with a quadratic error with respect to the original EPG 
without encoding or measurements. While our present 
construction addresses arbitrary linear decoherence, dif- 
ferent algebraic error structures may be tackled by mod- 
ifying the DD group. Notably, for dcphasing-dominated 



error processes, simpler DCGs based on 6 (vs 16) primi- 
tive gates suffices, which may be relevant to recently pro- 
posed fault-tolerant superconducting architectures under 
biased noise [20]. Our results provide the starting point 
for several generalizations. The restriction to driftless 
systems and the pulse shape assumptions may be relaxed, 
at the expenses of a more complex search for identifying 
distinct primitive gate sequences with the same leading 
error [1!)]. Systematic control faults can be further miti- 
gated by concatenating DCGs with composite pulses, at 
the expenses of longer control sequences. Conceptually, 
our analysis points to suggestive tradeoffs between avail- 
able error knowledge and error correctability in open-loop 
strategies. Ultimately, we believe that DCGs will further 
boost the practical significance of dynamical error control 
for quantum engineering and fault-tolerant computation. 

Support from the NSF under Grant No. PHY-0555417 
is gratefully acknowledged. 
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